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We investigate the dynamics of the FLRW flat cosmological models in which the vacuum 
energy varies with redshift. A particularly well motivated model of this type is the so- 
called quantum field vacuum, in which both kind of terms and constant appear in the 
effective dark energy density affecting the evolution of the main cosmological functions 
at the background and perturbation levels. Specifically, it turns out that the functional 
form of the quantum vacuum endows the vacuum energy of a mild dynamical evolution 
which could be observed nowadays and appears as dynamical dark energy. Interestingly, 
the low-energy behavior is very close to the usual ACDM model, but it is by no means 
identical. Finally, within the framework of the quantum field vacuum we generalize the 
large scale structure properties, namely growth of matter perturbations, cluster number 
counts and spherical collapse model. 
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1. Introduction 

The statistical analysis of various cosmological data (SNIa, Cosmic Microwave 
Background-CMB, Baryonic Acoustic Oscillations-BAOs, Hubble parameter mea¬ 
surements etc) strongly suggests that we live in a spatially flat universe that consists 
of ^ 4% baryonic matter, ~ 26% dark matter and ^ 70% some sort of dark energy 
(hereafter DE) which is necessary to explain the accelerated expansion of the uni¬ 
verse (see Refs. mm and references therein). Although there is a common agreement 
regarding the ingredients of the universe, there are different views concerning the 
possible physical mechanism which is responsible for the cosmic acceleration. 

The simplest dark energy candidate corresponds to a cosmological constant. In 
the standard concordance ACDM model, the overall cosmic fluid contains baryons, 
cold dark matter plus a vacuum energy (cosmological constant), that appears to fit 
accurately the current observational data and thus provides an excellent scenario to 
describe the observed universe. However, the concordance model suffers from, among 
other, two fundamental problems: (a) The fine tuning problem i.e., the fact that 
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the observed value of the vacuum energy density (pA = Ac^/SttG ~ 10“^^GeE^) 
is many orders of magnitude below the value found using quantum field theory 
(QFT),I^ and (b) the coincidence problem^ i.e., the fact that the matter energy 
density and the vacuum energy density are of the same order just prior to the 
present epoch, despite the fact that the former is a rapidly decreasing function of 
time while the latter is stationary. Such problems have inspired many authors to 
propose alternative candidates to the concordance ACDM, such as general dark 
energy models and modifications of the theory of gravity (for review see Ref. [5]). 

An alternative path that one can follow in order to alleviate the above problems 
is to consider a time varying vacuum, A(t). Within this framework we do not need to 
include new fields in nature nor to modify the theory of General Relativity. In this 
cosmological scenario the dark energy equation of state parameter w = Pde/pde, 
is strictly equal to -1, but the vacuum energy density is not a constant but rather it 
is a function of the cosmic time. Note that there is an extensive (olcPlSland nevPtSl) 
literature in which the time-evolving vacuum has been phenomenologically modeled 
as a function of time in various possible ways, in particular, as a function of the 
Hubble parameter. 

In this article we reconsider the so called quantum field vacuum (see next sec¬ 
tion). We would like to stress that this vacuum model was previously investigated 
by our team in many aspects. The layout of the article is the following. In section 2, 
we present the main cosmological ingredients of the quantum field vacuum model. 
In section 3 we provide the power spectrum and the linear growth of matter pertur¬ 
bations. In sections 4 and 5 we discuss the cluster number counts and the spherical 
collapse model. Finally, we draw our conclusions in section 6. 

2. Background expansion - Observational constraints 

The nature of the current vacuum modeP^^ is essentially connected with the renor¬ 
malization group (RG) in quantum field theory (QFT). In this context, the evolution 
of the vacuum is written as 

= + -Hi). ( 1 ) 

where Aq = A{Hq) = SHaLIq and e is interpreted in the RG framework as a 
“/3-function” of QFT in curved spacetime, which determines the running of the cos¬ 
mological constant. Regarding the dynamical role of e it has been found that since 
|i/| 1 (for a recent discussion see Ref. [13] and references therein) at low redshifts 

the model becomes almost indistinguishable from the cosmic concordance model. In 
other words the z/-parameter endows the vacuum energy of a mild dynamical evo¬ 
lution which could be observed nowadays and appears as dynamical dark energy. 
The low-energy behavior is thus very close to the concordance model, but it is by 
no means identical. 

Naturally, the next step here is to derive the Friedmann equations assuming of 
course a running vacuum. This procedure is perfectly allowed by the Gosmological 
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Principle embedded in the FLRW metric. In this context the Friedmann equations 
are written 


SirCptot = STrGpm + A = 3 





( 2 ) 


STrGptot = SttGp^ - a = -2H - (3) 

where the overdot denotes derivative with respect to cosmic time t. Notice that the 
Bianchi identities that insure the covariance of the theory, pose an energy exchange 
between vacuum and matter for G = const. 


Pm “t” d(l J- Pjyi — p\ . 


( 4 ) 


Combining equations ©, © and 0, we provide the basic differential equation 
that governs the dynamics of the Universe 

ij+^(l+cc„)iL^ = ^i:^A(iL). (5) 

Since we are in the matter era we set ujm = 0. 

Inserting eq. o into eq.© we can integrate the latter in order to obtain the 
Hubble parameter as a function of time: 


H{t) = Hq \ —- coth 

" 1 — V 


- Hq^ (Ha - ^)(1 - n) t 


( 6 ) 


where Ha = 1 — flm and Hq is the Hubble constan10. Using H = a/a the cosmic 
time, t{a), is given by 


t{a) 


3n 


1/2 


2 


sinh ^ 



where we have introduced 


n 


m 


l-v" 


Ha 


Ha — V 
l-v 


( 7 ) 

( 8 ) 


Inverting eq. © we easily determine the scale factor a = a{t). Substituting the 
cosmic time into ®, one can provide the normalized Hubble parameter 


E^{a) = = Ha + (9) 

where the cosmological parameters obey the standard cosmic sum rule, namely 
— 1 — H" 


^For the comoving distance and for the dark matter halo mass (see section 4) we use the tradi¬ 
tional parametrization Hq = lOOhkm/s/Mpc. Of course, when we treat the power spectrum shape 
parameter F we utilize h = hpianck = 0.673 [2]- 
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Concerning the evolution of matter density the situation is as follows. Combining 
eq. o, © and (HJ we arrive at pm + SHpm = Si/Hpm- Integrating the latter (using 
Pm = aHdpm/da) we find 

Pmia) = pmo , ( 10 ) 


where pmO is the matter density at the present time (a = 1), and therefore flm = 
Pmo/Pco, where pco = SHq/SttG is the current critical density. In fact, defining 
^m{a) = Praia)/Pc{a) it is easy to see, with the help of (fTUl) and the definition of 
E(a), that 


ia) 


E^ia) 


( 11 ) 


Lastly, upon inserting OT in 0 and integrating once more in the scale factor 
variable, we arrive at the explicit expression for the evolution of the vacuum energy 
density: 


A(a) = Ao + SttG 


^ PmO 
1-y 


- 3 ( 1 - 1 /) _ ^ 


( 12 ) 


Obviously, for ly = 0 the current time varying vacuum model reduces to the concor¬ 
dance A cosmology as it should. 

Using the above normalized Hubble parameter [see eq. ®]i a joint statistical 
analysis, involving the latest cosmological data [SNla,^ A(z) of BAOJl^and Planck 
CMB shift parameteiP^ is implemented. Since the explored cosmological models 
(Afic and ACDM) contain different number of free parameters, as a further sta¬ 
tistical test we use the {corrected) Akaike Information Criterion (AIC) relevant 
to our case {Ntot/k > 40),which is given, in the case of Gaussian errors, as 
follows: AIC = xlmin + ‘2-k, where Xt = XsNia + xIao + Xcmb is the overall 
function and k is the number of free parameters. From the statistical viewpoint, 
a smaller value of AIC points to a better model-data fit. Within this framework, 
we need to make clear that small differences in AIC are not necessarily signif¬ 
icant and therefore, it is important to derive the model pair difference, namely 
AAIC= AlCy — AICa;. The larger the value of [AAICj, the higher the evidence 
against the model with larger value of AIC, with a difference [AAICj > 2 indicating 
a positive such evidence and [AAICj > 6 indicating a strong such evidence, while 
a value < 2 indicates consistency among the two comparison models. For the con¬ 
cordance ACDM cosmology, if we impose 0 = 0 and minimizing with respect to 
rim we find rim = 0.291 ± 0.011 with xlmini^m)/dof ~ 567.5/586 (AICa^^ 569.5). 
For comparison, we provide the results of Planck: rim = 0.1426 ± 0.0025, with 
/ipianck = 0.673 ± 0.012 which implies rim = 0.315 ± 0.016. Notice that throughout 
the paper we set rim ,a = = 0.291. In the case of the Ajio model we find that 

the overall likelihood function peaks at rim = 0.282 ±0.012 and v = 0.0048 ±0.0032 
with Xt,mini^rn,v) ~ 563.8 (AICrg— 567.8) for 585 degrees of freedom.!^ It turns 
out that the current vacuum model appears to fit slightly better than the ACDM 
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the observational data. Still, the |AAIC| = |A1 Ca — A1C_rg| values (ie., < 2) indi¬ 
cate that the cosmological data are simultaneously consistent with the cosmological 
models. Note that for the rest of the paper we will restrict our present analysis to 
the choice of = (0.282,0.0048). 

Table 1. Statistical results from fitting SNIa+BAO+CMBs/^jj-^ data. The 
column indicates the model. The 2'^^ and 3^^^ columns provide the best fit 
parameters The and 5^^^ columns list the statistical significance 

of the fit. Finally, the 6^^ and 7^^ columns shows the estimated values of 
((78? <^c)- Notice that in section 4 we provide a relevant discussion regarding 6c. 

Model Qrn V Xt,min (78 (5c 

ACDM 0.291 ±0.011 0 567.5 569.5 0.829 1.675 

Krg 0.282 ±0.012 0.0048 ± 0.0032 563.8 567.8 0.758 1.644 


3. CDM Power Spectrum - growth factor 


The CDM power spectrum is given by P(fc) = {k), where T{k) is the CDM 

transfer function and n cii 0.9603 following the recent analysis of the Planck data.^^ 
Concerning the form of T{k), we use that of Eisenstein k HrP^ T{k) = 
where q = Lq = ln(2e ± l.Sq), e = 2.718 and Co = 14.2 ± j^§^- Here 
r is the shape parameteJi^ which is written as: P = Hm^pianck exp(—Hf, — 
\/2/ipianck Hb/ilm). The value of T, which is kept constant throughout the model 
fitting procedure, is estimated using the Planck results fib = Q with 
^Planck — 0.672. 

Another important quantity here is the rms fluctuations of the linear density 
field on mass scale Mh- 


a^{M,z) = al{z) 


k^‘+‘^T^{fl^,k)W‘^{kR)dk 

/“ k)W^{kRs)dk ’ 


crsiz) = crs 


Djz) 

D{0) ■ 


(13) 


where as = a{Rs,0) is the rms mass fluctuation on i?g = 8h~^ Mpc scales. At this 
point we would like to remind the reader that D{z) is the growth factor (see below), 
W{kR) = i{sinkR — kRcoskR)/{kR)^ and R — {8Mh /= flmPcO 
denotes the mean matter density of the universe at the present time {pmo = 2.78 x 
10^^Om/i^MQMpc“^). In this work we use the Planck prior, namely as ,a = 0.829. 
Of course, in order to use as properly along the current vacuum model we need to 
rescale the value of as by 


cs — <^8,A 


m 

Da{0) 


Po /“ fc"«+2T2(0„, k)W^{kRs)dk 
Po,A C fc"^+2r2(L!„,A, k)W^kRs)dk 


(14) 


where Pq/Po.a = [^m,A/^mf ■ Based on the aforementioned observational con¬ 
straints and eq. da we find as = 0.758 for the Arg model. In Table 1, one may see 
a more compact presentation of the cosmological parameters including the Sc which 
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is the linearly extrapolated density threshold above which structures collapse (see 
section 4). 

Now we focus on the basic equation which describes the evolution of matter 
fluctuations within the context of the previously described vacuum model. At sub¬ 
horizon scales the corresponding time evolution equation (in the linear regime) for 
the matter density contrast D = 5pmlPm, in a pressureless fluid, is given by (see 
Refs. [3 111]): 

b + {2H+ Q)t)-{A-KGpm-‘2HQ-Q^D = 0, (15) 


where in our case Q{t) = —K/pm- It becomes clear, that the interacting vacuum 
energy affects the growth factor via the function Q{t). In the case of non interacting 
DE models, [Q{t) = 0], the above equation (fT^ reduces to the usual time evolution 
equation for the mass density contrast Pqj. concordance A cosmology it is 
easy to prove that the solution of ea. dTSl) i^ 


DA{a) 


5nmE{a) r dx 

2 X x^E^{x) 


(16) 




Fig. 1. Right Panel: Comparison of the observed and theoretical evolution of the growth rate 
f{z)o-8{z). The dashed line corresponds to the ACDM (cts^a = 0.829) and the solid one is for A^q 
((jg = 0.758). Left Panel: The evolution of the growth index of matter perturbations. 


For the purpose of our study it is important to solve eq. m in order to investi¬ 
gate the matter fluctuation field of the vacuum model ([T]) in the linear regime. Here 
we present the main steps (for more details see Ref. [H]). First of all we change 
variables from t to a new one according to the transformation 


y = coth 


^ HoV (Ha v) t 


(17) 


Also from equations ®, © and m we get the following useful relations: 


y = 


1-v 
Ha — v 


E{a) , - 1 = 




Ha — V 


a 


- 3 ( 1 - 1 /) 


(18) 
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Utilizing (El) and m we find, after some non-trivial algebra, that equation ED 
becomes 


3 / 32 (y 2 _ 1 ) 2 ^ + k{l3)y{y^ - 1)D - 2[g{f3)y^ - i:{(3)]D = 0, 


a solution of which is 


D{y) = C{y^ - 1)-^ y-P" + 


1 3 1 


(19) 


( 20 ) 


^3/3 ' 2’ 2’ 3/3 2’ y^ - 1^ 

where j3 = \ — primes denote derivatives with respect to y, the quantity F 
is the hypergeometric function, fc(/3) = 2/3(6/3 — 5), y(/3) = (2 — /3)(3,9 — 2) and 
i/;(/3) = /3(4 — 3/3). Now inserting ea. dTSl) into ea. (|20)) and using eq.(| 8 l), we finally 
obtain the growth factor D{a) as a function of the scale factor: 


D{a) = Cv 


131 3 Ua 3 g 


( 21 ) 


where Ci oc C is an integration constant to be adjusted by an initial condition^ 

Furthermore, a crucial role in structure formation studies plays the growth rate 
of clustering which is defined as /(a) = ~ ^m(a) where 7 is the growth 

index. It has been shown that for those dark energy models which have a slow 
varying equation of state parameter, the growth index 7 is well approximated by 
7 ~ Refs. [Ml[53[55]), or 7 ~ 6/11 in the case where rc = —1 (ACDM 

model). Now since we know the analytical form of the growth factor D{a) one may 
directly compute the evolution of /(a). 

In order to test the performance of the Aug vacuum model at the perturbation 
level, we utilize the recent growth rate data for which their combination parameter 
of the growth rate of structure, f{z), and the redshift-dependent rms fluctuations 
of the linear density field, (Js{z), is available as a function of redshift, f{z)as{z). 
The total sample contains = 16 entries (as collected by Ref. [27] - see their Table 
1 and references therein). In the right panel of figure 1 we present the growth data 
(solid points) together with the predicted f(z)as(z) for the running vacuum model 
(solid line) and ACDM (dashed line). We find that the concordance A cosmology 
(with Dm = 0.291) reproduce the growth data with Xmin/^^f — 20.02/15 implying 
that the ACDM model can not simultaneously accommodate the Planck priors and 
the growth data (see also Ref. |1H])- On the other hand for the Afic model we have 
X^inldof ~ 8.42/15. Performing the AlC information criterion for small sample 
size, namely AlC = x^in + 2A: + we find AA1C= AICa — AIC^g — 11.6 

which means that the growth data favor the current running vacuum scenario. 

Let us finish this section with a relevant discussion regarding the growth index 
7 . As we have already described, we can express the linear growth rate of clustering 
in terms of the growth index. Specifically, from ~ ^m(a) we easily obtain: 

7 ( 2 ) ~ inQniiz)'' ^ where z = a~^{z) — 1. Based on the aforesaid observational 


^We normalize the growth factor at large redshifts {z ^ 1), namely D{z) ~ 1/(1 + z). 
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constraints we find that the growth index at the present time becomes 
which is somewhat higher with respect to the ACDM value, namely 7 ~ 6/11. In the 
left panel of figure 1 we show the growth index evolution for the h.RG and ACDM 
models, respectively. From the comparison we observe that the growth index of the 
A_rg vacuum model deviates from the concordance model. In particular, there is a 
visible deviation from above in all the redshift range. This deviation becomes at the 
level of 10 % for z < 1 . 




Fig. 2. Left Panel: Fractional difference 6J\f/Nf^ (solid line) in the cluster number counts between 
the A^rg model and the reference ACDM model. The dotted curves correspond to 20- Poisson 
uncertainties. Notice that we use the priors of Table 1. Right Panel: In the inner panel we provide 
the relative deviation (1 — A/Aa)% of the collapse factor for various vacuum models with respect 
to the A solution. In the outer panel we observe the evolution of the virial density. The lines 
corresponds to the following vacuum models: (i) Arqh (solid), (ii) Argci (dashed line, Us = 
0.002) and (iii) ArgC 2 (dotted line, Vs = —0.002). Notice that the solid points represent the 
concordance ACDM cosmology. 


4. Comparison with Cluster Halo Abundances 

Another choice in order to discriminate the models, that has been extensively used 
so far in the literature, is to estimate the theoretical predictions of the models for 
the cluster-size halo redshift distributions and to confront them with the data. The 
halo abundances predicted by a large variety of time varying vacuum models have 
been compared with those corresponding to the ACDM model.EiHinillll yvb utilize 
the Press and Schechtei^ formalism, which determines the fraction of matter that 
has formed bounded structures as a function of redshift. Mathematical details of 
such a treatment can be found also in Ref. |13) ; here we only present the basic ideas. 
The number density of halos, n(M, z), with masses within the range (M, M -b JM) 
is given by: 

p dlntj”^ 


n(M, z)dM 


M dM 


f{a)dM , 


( 22 ) 
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where p = pmo is the mean background mass density (see section 3). In the original 
Press-Schechter (PSc) formalism, f{a) = /psc(o’) = y^2/TT{Scjcr) exp{—S‘^l2a^), 
5c is the linearly extrapolated density threshold above which structures collapse,!^ 
while z) is the mass variance of the smoothed linear density field, extrapolated 

to redshift z at which the halos are identified [see ea. dTJt ]. The mass variance 
depends on the power-spectrum of density perturbations in Fourier space, P{k), 
for which we use here the CDM form according to Ref. EQ], and the values of 
the baryon density parameter, the spectral slope and Hubble constant according 
to the recent Planck results.^ Although the original Press-Schechter mass-function, 
/pSc, was shown to provide a good first approximation to that provided by numerical 
simulations, it was later found to over-predict/under-predict the number of low/high 
mass halos at the present epoch.l^l^ More recently, a large number of works have 
provided better fitting functions of /(cr), some of them based on a phenomenological 
approach. In the present study, we adopt the one proposed by Reed et al. [34]. 

In order to compare the mass function predictions of the different cosmological 
models, it is important to use for each model the corresponding value of erg and 6c 
(see Table 1). It is well known that for the usual A cosmology we have 6c — 1.675 
(see Ref. [35]). For the Arc model it has been found that 6c — 1.644 (see Ref. [I3]). 
Given the halo mass function from eq. (1221) we can now derive an observable quantity 
which is the redshift distribution of clusters, Af{z), within some determined mass 
range, say Mi < M/h~^MQ < M 2 = 10^®. This can be estimated by integrating 
over mass the expected differential halo mass function, n{M, z): 

M{z) = — / n{M,z)dM, —=n,r{z) — , r{z) = c tTTT > (23) 

dz Jmi dz dz Jo H{u) 

where dV/dz is the comoving volume element and fig is the solid angle. 

In left panel of figure 2, we present the fractional difference between the A/jg 
model and ACDM namely, 6J^/N\, with masses in the range h~^ < M/Mq < 
10^^ h~^. Recall that we denote the deviations of the number counts of a given 
vacuum model with respect to the ACDM as JAC = A/"—A/a. The dotted lines shown 
in the left panel of figure 2 correspond to 2a Poisson uncertainties, which however 
do not include cosmic variance and possible observational systematic uncertainties, 
that would further increase the relevant variance. Specifically, the prediction for 
5fJ'/J\f is negative for all points, even for those in the ±lcr band. Thus, the running 
vacuum model tends to predict a smaller number of counts as compared to the 
ACDM. The fractional decrease can be as significant as 30 — 60%. We also find that 
the redshift variation of the differences between the A/jg cosmology and ACDM 
model is mainly affected by variations in the values of erg and Jc (for more discusion 
see Ref. [E]). Therefore, we have verified that there are observational signatures 
that can be used to differentiate the A^g model from the ACDM and possibly from 
a large class of DE models. 
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5. The spherical collapse model 

In this section we generalize the spherical collapse model within the variable A{H) 
cosmological model, in order to understand non-linear structure formation in such 
scenarios and investigate the differences with the corresponding expectations of the 
concordance ACDM cosmology. Practically, one may start from the Raychaudhuri 
equation which is valid either for the entire universe or for homogeneous spherical 
perturbations [by replacing the scale factor with radius R(t)] 


R 47rG/ N 

^ ^ {Pms ~ ^PAs) , (24) 

where pms and pAs refer to the corresponding values of the matter and vacuum 
energy densities in the spherical patch susceptible of ulterior collapse. 

Now let us first define some basic quantities of the problem. In particular, we 
call at the scale factor of the universe where the overdensity reaches at its max¬ 
imum expansion (i.e. when R = 0) and ac the scale factor in which the sphere 
virializes, implying that a cosmic structure has formed. Also, Rt and Rc denote the 
corresponding radii of the spherical overdensity, the former being the turnaround 
(or “top hat”) value at the point of maximum size, and the latter refers to the 
eventual situation when the sphere has already collapsed and virialized. We would 
like to remind the reader that due to the coupling between the running vacuum 
and matter one would expect that the matter density in the spherical overdensity 
should obey the same power law as the background matter Pm{a) oc (see 

eqUni). Thus, Pms oc R denotes the matter density in the spherical patch. 

Regarding, the vacuum energy density in the spherical region, pAs we focus on two 
situations: (i) the vacuum energy remains homogeneous and only the corresponding 
matter virializes (this scenario holds for the A cosmology); and (ii) the case with 
clustered vacuum energy, assuming that the whole system virializes (both matter 
and vacuum components). 

In the regime where the vacuum is allowed to cluster we have pas{R) = 
As(i?)/87rG. Therefore, in such a situation it could be possible, on non-linear scales, 
to have an interaction between dark matter and dark energy with a different jz than 
the background value. We also assume that the general functional form that de¬ 
scribes the behavior of the vacuum energy density inside the spherical perturbation 
obeys a similar equation as that of eq. (HU): 


A.(i?) 


Ao -f SttG 


Pms,t 
1 - Vs 



(25) 


where Vg is not necessarily equal to the background v. 

Using the basic differential equations [see eq. m and (f24|l ] and performing the 
following transformations x = ajat and y = R/Rt we arrive at 


X 


2 


H^nm,t + rx^Hx)] 


(26) 
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and 


y 

y 



C _ g PAs 

y Hm,t 


(27) 


where Notice that Um.t = is the matter density param¬ 

eter at the turnaround epoch (see eqE]). 

It is important to point that in order to obtain the above set of equations we 
have used the following relations: 


and 


Pms 


I{x) 


PA 1 - 1 - h'roap^^~''^ 
PKt 1 -I- i/roGf 


PA,t 
r = - 

Pm,t 


3(l-y) , 

+1 - 

ilm I -I' 



3{l-u) 

t 


(28) 

(29) 


(30) 


where tq = flA/^m, Pm,t, PA,t are the matter and the vacuum energy density at 
the turnaround epoch which satisfies Ha,* = 1 — (for definition see ea fTl]) . 
Of course we need to say that in order to derive eqs.([29|) and (1^ we have used 
eqs.([8|) and uni). The matter density in the spherical region at the turn around 
time is given with respect to the background matter density at the same epoch 
Pm,t as Pms,t = CPm,t- The parameter C, is the density contrast at the turnaround 
point. It is interesting to mention that in the case of the Einstein-de Sitter model 
(Om = = 1 and 1 ^ = 0) the solution of the system formed by ea. (l26)l and eq.([2ll) 

reduces to the well known value of the density contrast at the turnaround point: 
^ = (^)^, as it should. 

For bound perturbations which do not expand forever, the time needed to re- 

(7]), we obtain 

~ 2 sinh“^ 


collapse is twice the turn-around time, tc — 2t(. Therefore, from eq. 
the relation between the ac and at sinh“^ 


ro a, 


3(1-1.) 


3(1-!.) 

roOt' 


which is well approximated by Zt — 1.523zc-1-0.8. Therefore, considering that clusters 
have virialized at the present epoch, Zc — 0, the turnaround redshift is zt — 0.8. 
If we assume that galaxy clusters have formed close to the epoch of z^ ^ 1.5 then 
we obtain zt ^ 3. To this end we verify that the ratio between the scale factors 
converges to the Einstein de Sitter value {ac/at)ARG — (1 + Zt)/{1 + Zc) = 2^/^ at 
high redshifts owing to the fact that the matter component dominates the Hubble 
expansion. 

In order to provide the virial theorem of the Arg model we need to generalize the 
Layzer-Irvine equation, which describes the flow to virialization.l^ If we take into 
account that the matter is exchanging energy with the vacuum then the modified 
virial theorem become^^ 


(2 - 3i/)r -h (1 - 6i^){Ug - 2Ua) = 0 


(31) 
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where Ug = 3GM'^/5R. The vacuum potential energy Ua is written as 


A{a)MR^ 
10 


Ua — { AoMR^ I iTrGvsMp ^R^ _ iTTGvsMp„^s,tR-^^'^''^ 


10 


Homogeneous 

Clustered 


(32) 


where M is the mass inside the spherical overdensity. Utilizing the observational 
constraint v = 0.0048, the deviation from the usual virial condition is ^ 2 — 3%. In 
this context, combining the virial theorem and the energy conservation (T'c + 17g,c + 
Ua,c = UG,t + UA,t) at the collapse and at the turn around epochs we reach to the 
following condition: 

qi{v)UG,c + q 2 {v)UA,c = UG,t + UA,t, qi{v) ^ q 2 {v) = . (33) 

2 — 6u 2 — 6v 

Clearly, for v = 0 the above equations boil down to those of the A model. 


5.0.1. Homogeneous vacuum 

In the case of homogeneous vacuum (hereafter A ugh) model we have pAs{ct) = 
PA (a) = A(a)/87rG. Within this framework, inserting ea. (l2^ into ea. (l27l) . we obtain 


y = 


H?nra,t 


ly 


S^-^ryI{x) 


(34) 


The solution for ^ is provided only numerically by integrating the main system of 
differential equations, leasIMland IMll . using the boundary conditions: (dy/dx) = 0 
and y = 1 at a: = 1. However, we find that a reasonably accurate fitting formula for 
C is given by 

2 








(a;i,a;2,a;3) = (0.79,0.26,0.06) (35) 


where w{a) = —1 —Using egs. (15^ . (1551) and Ug = 3GM‘^/5R we can obtain a 
cubic equation that relates the ratio between the virial Rc and the turn-around outer 
radius Rt the so called collapse factor (A = Rc/Rt), q 2 {i^)ncX^ — {2 + nt)X-\-2qi{i') = 
0 , where 


«c.t = 


A(ae,t) 

47rGpm,tC 


= «o + 


C(i-^) 


-3(1-!.) 


- 1 


(36) 


with no = ——. The viable solution (0 < A < 1) of the above cubic equation 
is 

'-27r' 


, 2dU3 
A =--—cos 


(~ 3 ~) ’ d = \/+ xl, 6 = cos ^{xi/d) ( 37 ) 

where xi = —21qi, X 2 = — and = 4 ■ Of course in the case 

of z. = 0 the above expressions get the usual form for A cosmologj^^^ while for an 
Einstein-de Sitter model (17^ = 1) we have A = 1/2. 
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We find that the collapse factor lies in the range A ~ 0.48 — 0.50 in agreement 
with previous studies|-]^g inner panel of figure 2 we plot the relative 
deviation (solid line) of the collapse factors Xargh i^c) for the current vacuum model 
with respect to the A solution Aa(zc). Obviously, the deviation from the ACDM case 
is small ~ —2%. In the outer panel of figure 2 we present the evolution of the density 
contrast at virialization (solid curve) 


_ Pms,c C / 1 

Pm,c A^ \(lt J 


(38) 


We verify, that the density contrast decreases with the virialization redshift z = Zc 
and at very large redshifts it tends to the Einstein-de Sitter value (A„ir ~ IStt^), 
since the matter component dominates the cosmic fluid. Following the notations of 
Ref. [35j , we obtain an accurate fitting formula to Ayir (within a physical range of 
cosmological parameters and for z < 2) 


A„ir(a) — IStt^ [l + e ©'’(a)] , 0(a) = O^^(o) - 1, (39) 

where e = 0.40 — 2S.2v + 500j^^ and b = 0.94 + S9.6v. From the right panel of figure 
2 we observe that the virial density Ayir of the Augh model is somewhat lower 
with respect to that of ACDM model, namely the relative difference can reach up 
to ~ —4%. As an example, assuming that clusters have formed prior to the epoch of 
Zc ~ 1.5 (zt 3) we find (C, A„ir)ARGH - (5-66,183.9) and (C, A„ic)A - (5.65,190). 
To conclude this discussion we would like to stress that in the homogeneous case 
the dark energy (in our case vacuum) component flows progressively out of the 
overdensitjSH^ and hence energy conservation cannot be applied. Such violation 
is however only important for large values of |iz| ~ 0(10“^) and at very late times, 
when the vacuum energy dominates the cosmic fluid. 


5.0.2. Clustered vacuum 


In this section we assume that the vacuum clusters along with the dark matter. As 
we have already mentioned the vacuum inside the spherical patch is given by ea. (l25p . 
In the current study we restrict our analysis to = 0.002 (hereafter Argci'- dashed 
line) and Vs = —0.002 (hereafter AjiGC 2 - dotted line). Substituting y = R/Rt into 
eq. dSSl) we have 


As(y) 


Aq + SttG 


CPm,t 
I - Vs 


-Hi-u,) _ ^ 


(40) 


where the vacuum energy density is pAsiv) = As(y)/87rG. The merit of the latter 
assumption is that it allows an analytical solution to the system of eas. (l26)) and 
(HZl). Indeed inserting the above form of PAsiu) into eq (071) we obtain 


y = - 




2 


(1 - _o(r- 

(l_^^)y2-3.. 



(41) 
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and upon integration we find 


y^=H^nm,t[Piy,C)+C], (42) 

where C is the integration constant and P{yX) = [i-i, 

Notice that the boundary conditions, (dy/dx) = 0 and y = 1 at a: = 1, imply that 
C = —P{1X)- Combining ea. d^ with the background equation ([26|) the solution 
of the system is 


dy 


dx 


/o v'^(y,C)-^(l,C) 




-l+3i^ 


+ rx^I{x) 


(43) 


Furthermore, if the vacuum energy participates in the virialization then the poten¬ 
tial energy of the vacuum is given by the second branch of ea. (l32|) . In this case, 
using simultaneously eos. d^ . d^ and Uq = 3GM‘^/5R the collapse factor obeys 

q2Xs)[no - /(j^s)]A^ - A{no, Vs)X + g{vs)X^''‘ + 2qi{vs) = 0 , (44) 


where 

= 


2k, 


gips) = 


lQVsq2Xs) 


A{no,Vs) = 2 +no - f{vs) + 


gXs) 


l-uX {l-ys){2-iysy ' qXnsY 

Finally, solving eas. (l4^ and (lUll . we can estimate the density contrast at virial¬ 
ization from ea. (l38|) . Similarly, as in section 6.2.1, the fitting formulas for C as well 
as for /Xyir are given by ea. d^ and ea. dMl) . The corresponding parameters of the 
approximated formulas are 


(a;i,W2,W3) 


(0.62, -0.08,0.06) Q<vs< 0.003 
(0.82,0.22,0.06) - 0.003 < < 0 


(45) 


(h (0.94-h90j^«, 0.40-46 i/,-H 500j^2) - 0.002 < < 0.003 , . 

^ t (0.94-h95jz„ 0.31-129i^s-f 500i/2) - 0.003 < <-0.002 ^ ’ 

In this framework the collapse factor obeys 0.46 < A < 0.52. Also from figure 
2 (see inner and outer panels) we see that the largest positive deviation of the 
collapse factor occurs for the Argc 2 model (dotted line). This implies that ArgC 2 
model forms more bound systems than the concordance ACDM model (solid points) 
and thus the corresponding cosmic structures should be located in larger den¬ 
sity environments. Indeed, at the cluster formation epoch Zc — 1.5 we obtain 
((, Ayir)ARGC 2 — (5.66,202). The opposite situation holds for the Arggi model. 
From figure 2 it becomes obvious that the size and of the cosmic structures 
which are produced in the Argci model (dashed line) are remarkably close to that 
predicted by the Argh vacuum cosmology, and therefore the impact of the vacuum 
energy on the spherical collapse is very small in the clustered case as long as iZg is 
positive. In other words we find that both Argh and Argci models are equivalent 
at the background and perturbation levels. Finally, at the epoch of Zc — 1.5 we find 
(C, A„ir)AHGci — (5.70,184.2). For comparison we provide the density pair of the 
Argh model, namely (C, A„ir)AflGH - (5.66,183.9). 
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6. Conclusions 

In this review article we have studied the overall dynamics of the FLRW flat cos¬ 
mological models in which the vacuum energy varies with the Hubble parameter, 
namely A{H) = Ag -I- — Hq). First we have performed a joint likelihood anal¬ 

ysis in order to put constraints on the main cosmological parameters by using the 
current observational data (SNIa, BAOs and CMB shift parameter together with 
the growth rate of galaxy clustering). We have shown that the A{H) model fits 
slightly better the observational data than that of the traditional A cosmology. In 
particular, we have found that the ACDM model can not simultaneously accommo¬ 
date the Planck priors and the growth data implying that this kind of data favor the 
A(iL) vacuum scenario. Subsequently we have investigated the nonlinear regime and 
considered the predicted redshift distribution of cluster-size collapsed structures as 
a powerful method to distinguish the A{H) and ACDM cosmological scenarios. Fi¬ 
nally, we have generalized the properties (virial theorem, collapse factor, virial and 
turnaround densities) of the spherical collapse model in the case when the vacuum 
energy is a running function of the Hubble rate, A = A(H). Overall, we have found 
that the virial density contrast is affected by the considered status of the vacuum 
energy model (homogeneous or clustered). 

Acknowledgments: I would like to thank J. Sola, M. Plionis, J. A. S. Lima 
D. Polarski, N. E. Mavromatos, L. Perivolaropoulos and A. Gomez-Valent for the 
recent collaboration in some of the work presented in this review article. 
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